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ABSTRACT 

Context. Most of the observed emission lines and continuum excess from young accreting low mass stars (Classical T Tauri stars - 
CTTSs) take place in the star-disk or inner disk region. These regions have a complex emission topology still largely unknown. 
Aims. In this paper the magnetospheric accretion and inner wind contributions to the observed permitted He and H near infrared (NIR) 
lines of the bright southern CTTS RU Lupi are investigated for the first time. 

Methods. Previous optical observations of RU Lupi showed a large Ho- profile, due to the emission from a wind in the line wings, 
and a micro-jet detected in forbidden lines. We extend this analysis to NIR lines through seeing-limited high spectral resolution 
spectra taken with VLT/ISAAC, and adaptive optics (AO) aided narrow-band imaging and low spectral resolution spectroscopy with 
VLT/NACO. Using spectro-astrometric analysis we investigate the presence of extended emission down to very low spatial scales (a 
few AU). 

Results. The Hel /i 10830 line presents a P Cygni profile whose absorption feature indicates the presence of an inner stellar wind. 
Moreover the spectro-astrometric analysis evidences the presence of an extended emission superimposed to the absorption feature 
and likely coming from the micro-jet detected in the optical. On the contrary, the origin of the Hydrogen Paschen and Brackett 
lines is difficult to address. We tried tentatively to explain the observed line profiles and flux ratios with both accretion and wind 
models showing the limits of both approaches. The lack of spectro-astrometric signal indicates that the HI emission is either compact 
or symmetric. Our analysis confirms the sensitivity of the Hel line to the presence of faint extended emission regions in the close 
proximity of the star. 

Key words, stars: individual: RU Lupi - stars: winds, outflows - accretion, accretion disks - line: formation - techniques: high 
angular resolution 



1. Introduction 

Classical T Tauri stars (CTTSs) are young solar and lower mass 
stars presenting a rich emission line spectrum in excess of the 
continuum from the stellar photosphere. This emission includes 
Hydrogen lines, low energy forbidden lines and Helium lines 
together with a continuum excess - the veiling. 

Optically and NIR forbidden emission lines, such as the 
[Oi], [S n] and [Fe ii] lines, have emission regions large enough 
to be r esolved and are typically associated with outflowing 
jets (e .g. lBacciotti et al.Ll2000HDougados et alLlIioOOUPvo et all 
I2OO6 I). On the other hand, the atomic permitted lines, such as 
the Hydrogen lines, are more difficult to interpret. The energy 
released by the accretion of matter onto the star can power 
the permitted line emission that is believed to arise from mag- 
netospheric accretion columns and/or accretion shocks (e.g. 
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* Based on observations collected at the European Southern 
Observatory, Chile (ESO Programme 71.C-0703). 



iMuzerolle et"an 120011: iBouvier et all l2007l) . The fact that the 
emission component of permitted atomic lines is attributed to 
infalling gas is supported by the presence of inverse P Cygni 
profiles, and by the correlation of their luminos ities with the ac- 
cretion luminosity, derived from the veiling ( MuzeroUe et al.L 
ll998blH) . On the other hand, the blueshifted absorption fea- 
tures, when pres ent, are an indicator of high veloc ity inner 
winds (e.g. Edwa rds et all Il987t iNaiita et aU |2000|) . Indeed, 
several recent spectroscopic studies suggest that the emis- 
sion of the permitted HI and He lines can actually include 
a contribution from b oth inflowing and outflowing gas (e.g . 
'Beristain et al.', '2001'; Folha & Emerson, '200 li iTakami et all 
2001; Whelan et al., 2004; Edwards et al, 2006). This has moti- 
vated hybrid models where the Hydrogen e mission is due t o both 
the magnetospheric ac cretion and a wind (lAlencar et all 120051: 
iKurosawa et al.Ll2006l) . 

The fact that the collimated, spatially resolved jets, traced 
at large distance from the source by optical, near-infrared 
forbidden lines, are generated in circumstellar regions (less 
than 10 AU from the source) is predicted by theoretical models 
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dAnderson et al.Ll2003tlFerreira et al.Ll2006h and observationally 
proved by the frequent detection of blueshifted absorption 
featur es in a number of s trong permitted lines (Ha, Na, D, Call, 
Mgll, iNaiita etal] |2000|) . formed close to the star Signatures 
of outflowing gas revealed by absorption featur es can also be 



searched in the profile of the Hel /110830 line jTakami et al 



12002- Edwar ds et all llOOl iDupree et"ail l2005h . Edwards et alj 
(2006) analysed a sample of 3 1 T Tauri stars and showed that the 
Hel line is much more sensitive than the usually used Ha line 
to trace inner winds in absorption (in 71% of the sources was 
detected an absorption featur e in the Hel fine again st only 10% 
of detections in the Ha line). 'Be ristain et al] (1200 ll) showed that 
the diagnostic power of permitted Helium lines in tracing inner 
winds, lies in their high excitation potential (20-50 eV), that 
restricts the line formation to a region of either high temperature 
or close proximity to a source of ionising radiation. 

In order to investigate the reliability of such scenarios, we 
have observed the circumstellar region of the T Tauri star RU 
Lupi in the NIR wavelength range. 

This so urce is one of the most active and well known T 
Tauri stars dLamzinetali 119961: IStempels & Piskunovi l2002t 
iHerczeg et aU 2005 1) andit is a good target to investigate the in- 
ner winds and accretion flows down to very sm all spatial scales . 
RU L upi is located at a distance of only 140 pc dde Zeeuw et al.L 
1199 9). allowing us to reach small spatial scales with high an- 
gular resolution instrumentation and/or specialised techniques 
of data analysis. In particular, NACO at the VLT was used to 
take high angular resolution images and low spectral resolu- 
tion spectra. With VLT/ISAAC seeing limited and high spectral 
resolution spectra were obtained. A spectro-astrometric analy- 
sis was applied. Spectro-astrometry is a powerful technique to 
obtain positional information on the region originating the emis- 
sion in the components o f a line profile even in seeing-limited 
conditions. Taka mi et alj (|2001), for example, using spectro- 
astrometry, found that the blueshifted and redshifted wings of 
the Ha line emitted by RU Lupi show a spatially symmetric off"- 
sets with respect to the source, in the same direction of the jet 
trace d, on higher spatial sc ales, by the forbidden [O i] and [S n] 
lines. IWhelan et all (|2004|) detected a similar extended emission 
in the wings of the Pa/S line in three T Tauri sources (DG Tau, 
V536 Aql, andLkHa 321). 

In this paper we present the results obtained from such a 
study and, in particular, a detailed analysis about the origin of 
the Paschen and Brackett HI lines and the Hel /110830 line. 
The paper is structured as follows: in Sect. |2] we present the 
different set of observations and the data reduction process: in 
Sect.[3]the observed line profiles and fluxes are analysed: more- 
over we present the concepts underlying the spectro-astrometry 
technique and the obtained position spectra: finally in Sect. |4] we 
discuss the origin of permitted H and He lines by means of a 
model for wind and accretion, and of the results obtained with 
spectro-astrometry. In the Conclusions (Sect. |5]l we summarize 
our findings. 



2. Observations & Data Reduction 

Observations of RU Lupi were carried out in between April and 
June 2003 at the ESO Very Large Telescope. The first set of 
spectra was acquired at the UT4-YEPUN telescope using the in- 
frared camera and spectrometer CONICA, working in cascade 
to the adaptive optics (AO) system NAOS, and the 86 mas width 
slit. The AO system allows one to obtain a high spatial reso- 
lution of ~a.'08-0'.'12 (full width at half maximum - FWHM 



- of a point source), sampled by a pixel scale of 0'.'054. The 
second set of data was taken at the UTl-ANTU telescope us- 
ing ISAAC with the Of! 3 slit to obtain a high spectral reso- 
lution (R~10000: spectral sampling: ~14 kms"': pixel scale: 
U.' 147/pixel), while the spatial resolution is limited by the seeing 
(FWHM~0'.'6-0'.'8). The NACO observations cover all the J and 
H bands wavelength ranges (R~400 and R~1500 respectively) 
while the ISAAC spectra centered at 1.08, 1.27 and 1.588 //m 
cover very short wavelength segments. Both the NACO and the 
ISAAC spectra were obtained at two slit position angles (217° 
and 127°), i.e. parallel and perpen dicular to the jet direction as 
determined by iTakami et al.l (1200 1). 

In addition two sets of images of the source in the 
narrow band filters centered around the HeI/110830 and 
[FeII]/tl.25,1.64yL(m lines, and around their near continuum (at 
1.09, 1.24 and 1.75 //m), were taken with NAOS/CONICA. 

The information about the spectral observations and the 
images are summarized in TablT] 

The spectra were reduced using standard IRAF tasks and 
Eclipse routines from the ESO pipeline. The images were cor- 
rected for bad pixels and cosmic rays and for ghosts, then flat- 
fielded, and sky-subtracted. Wavelength calibration was per- 
formed using the spectra of the arc lamps for the NACO data 
and the sky OH lines for the ISAAC spectra. The systemic mo- 
tion of the object was subtracted using the r adial velocity in the 
local standard of rest frame, determined bv ITakami et al.l(l200lb 
through the Li /16707.815 absorption line on the stellar atmo- 
sphere, Vlsr=8+2 kms"'. The spectra were corrected for the at- 
mospheric transmission and for the efficiency of the camera, and 
then flux calibrated using a telluric spectro-photometric stan- 
dard. The typical uncertainty for the flux calibration is of ~5%. 
Since the observed telluric stars are of type B, the absorption HI 
lines were subtracted before the correction. Unfortunately the 
Brio, Brl3 and Brl5 absorption lines were deeply blended with 
the atmospheric features in the telluric spectra. In the NACO 
spectra, because of the low resolution, it was not possible to 
perform a de-blending and subtract these HI absorption lines 
in the spectrum of the telluric. Thus the fluxes and the equiv- 
alent widths of the BrlO, Brl3 and Brl5 lines measured from 
the NACO spectra are only upper limits. In the ISAAC telluric 
spectra, on the contrary, the HI lines were de-blended from the 
atmospheric features and subtracted. Even at the higher spec- 
tral resolution of the ISAAC data, however, this process is not 
very accurate and implies an oversubtraction. For this reason the 
fluxes, the FWHM and the equivalent widths of the BrlO, Brl3 
and Brl5 lines inferred from the ISAAC spectra should be con- 
sidered as lower limits. 

The NACO images were reduced using standard procedures: 
removal of bad pixel and cosmic rays, dark subtraction, flat- 
fielding and sky-subtraction. Then the images acquired with the 
narrow band filters centered on the continuum emission were 
subtracted from the images filtered on the line emission (Hel and 
[Fe n] ). The subtraction was performed by: (i) choosing images 
with similar PSF, to minimize the effects of a different distribu- 
tion of the emission in the different filters: (ii) normalizing the 
images for the source flux estimated trough a gaussian bidimen- 
sional fit. 



3. Results 

Analysis of the NACO imaging data showed that no physical 
companion is detected in the continuum frames in the limit of 
our spatial resolution (~ 14 AU). No extended emission structure 
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Table 1. Instrumental setup during observations 



Instrument 


"/pix 


Slit 


Date 


A ipim) 


R 


PA 


T„„xN,, 


NACO 


0.013 


n.a. 


20/25-5-03 


NB 1.08, 1.09 


n.a. 


n.a. 


10.8s X 12 


imaging 








NB 1.24, 1.26 
NB 1.64, 1.75 






9sx 12 
9sx 12 


NACO 


0.054 


0'.'086x40" 


16-6-2003 


I: 0.91-1.37 


400 


217°, 127° 


18sx 10 


spectroscopy 






22-4-2003 


H: 1.50-1.84 


1500 






ISAAC 


0.147 


0'.'3x 120" 


26-6-2003 


SZ: 1.06-1.1 
J: 1.24-1.3 
H: 1.54-1.62 


11500 
10500 
10000 


217°, 127° 
217°, 127° 
217° 


300s X 4 
300s X 4 
250s X 2 



is detected in continuum subtracted emission line frames at Hel 
and [Feii]/ll.25,1.64yum. This constrains the presence of wind 
emission in the Hel and [Fe n] lines over a region smaller than 
14 AU. 

The NACO low resolution spectra show strong hydrogen 
recombination lines in emission (BrlO to Brl9 and Pa/3, Pay) 
and the Hel 1.0830//m line (see Figffland Tab.©. The ISAAC 
medium spectral resolution spectra allow us to study in more de- 
tail the spectral profiles of the Hel, Pa/?, Pay and Br 13, 14, 15, 
16 lines (see Fig|2]and Tab. [3]l, as well as their spatial displace- 
ment with respect to the continuum emission through spectro- 
astrometry (see Sect. 13.31 ). No [Feu] lines are detected at 1.25, 
1.32, 1.53, 1.60 and 1.64 fim both in the NACO and in the 
ISAAC spectra. The lack of [Fe ii] emission is not expected in 
this source. In fact, studies of HH jets in the optical and in the 
NIR range dPodio et"aLll2006 l) showed that optical jets emit also 
in the NIR [Fe n] lines, unless the density in the jet beam is very 
low (n^ < 10^ cm"^). In RU Lupi, indeed, on one hand strong and 
spatially extended emission in the optical forbidden line s ([S ii], 
[Oi], and [Nil]) was observed by iTakami et al.l (1200 Ih : on the 
other hand, electron densities higher than 10^-10^ cm"^ were es- 
timated by the same authors based on the [S ii] lines equivalent 
width ratio. Therefore, iron lines should be present in the spec- 
tra. Possible reasons for the non detection of Fe emission could 
be the fact that a short integration time was used, or that a strong 
depletion of iron caused an emission fainter than the predicted 
one. In fact, an increase of depletion near t o the source has bee n 
observed in other jet sources as, e.g. HH 1 dNisiniet aI.Ll2005h . 

Table 2. Equivalent widths (EW) and dereddened line fluxes (F) 
from NACO spectroscopy. 



A 


Line 


EW 


F 


fim 




km s"' 


10"" erg s"' cm"^ 


1.0830 


Hel 


452 


40 ±5 


1.0941 


Pay 


457 


40 ±5 


1.2822 


Pa/3 


434 


32 ±3 


1.5265 


Br 19 


9 


0.5 ± 0.2 


1.5346 


Br 18 


14 


0.8 + 0.2 


1.5443 


Br 17 


20 


1.1 +0.2 


1.5561 


Br 16 


19 


1.1 ±0.1 


1.5705 


Br 15 


46 


2.7 ± 0.3 * 


1.5885 


Br 14 


31 


1.8 ±0.2 


1.6114 


Br 13 


55 


3.1 ±0.3* 


1.6412 


Br 12 


48 


2.7 ± 0.3 


1.6811 


Br 11 


63 


3.5 ±0.3 


1.7367 


Br 10 


78 


3.9 ± 0.3 * 



* The equivalent widths (EW) and the fluxes (F) of the BrlO, Brl3 and 
Br 15 line are only upper limits to the real flux, due to the presence of HI 
lines blended with atmospheric features in the spectrum of the telluric 
standard (see Sect.|2]l. 
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Fig. 2. ISAAC spectra of the Hel /110830, Pay and Pa^ lines. 
The lines are normalized to the continuum intensity. 



Table 3. Line properties from ISAAC spectroscopy. 



A 


Line 


^ peak 


AV 


EW 


F" 


fjm 




km s"' 


kms"' 


km s"' 


10"" ergs"' cm 


"1.0830 


Helab 


-200 




133 




"1.0830 


Hel em 


70 




581 


47 ± 6 


1.0941 


Pay 


-5 


206 


475 


40 ±5 


1.2822 


Pa^ 


-11 


214 


562 


38 ±3 


1.5561 


Br 16 


-31 


179 


23 


1.2 ± 0.2 


* 1.5705 


Br 15 


-16 


128 


20 


1.0 ±0.2 


1.5885 


Br 14 


-22 


185 


43 


2.2 ± 0.2 


*1.6114 


Br 13 


-19 


159 


44 


2.7 ± 0.3 



" The fluxes (F) are corrected for dust extinction (see Sect l3.lt . 

Parameters relative to the absorption (ab) and emission (em) part of 
the Hel line profile (see Fig.|2j- 

* Due to the presence of HI lines blended with atmospheric features 
in the spectrum of the telluric standard, the parameters inferred for the 
Brl3 and Brl5 lines are only rough estimates (see Sect.|2ll. In particular 
the width, AV, and the flux, F, are lower limits. 



3.1. HI emission lines. 

The availability of complementary data sets (see Tab. [U allows 
us to investigate the origin of HI lines through both the mea- 
sured line fluxes and the line profiles. The characteristics and 
the parameters of the lines profiles can be retrieved from the 
ISAAC spectra thanks to the high spectral resolution. The cen- 
tre (Vpeak), the FWHM (AV), the equivalent width (EW), and 
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Fig. 1. Low resolution NACO spectrum of RU Lupi in the J and H bands. The Hel /110830 line and the HI emission lines from 
the Paschen and the Brackett series are showed. The dotted Unes indicate the wavelength position of the [Fe ii] lines at 1 .25 and 
1 .64yum, that were not detected in our spectra. 



the dereddened flux (F) of the HI lines which are in the spec- 
tral range covered by the ISAAC spectra (Pa^, Pay and Br 
13, 14, 15 and 1 6) are summarized in Tab. |3] The classifica- 
tion proposed by iReipurth et al.l d 19961) for th e Ha Une profiles 
was al ready used for N IR Hydrogen lines by Folha & Emerson 
( 1200 lb and IWhelan et al. ( 2004 ). According to this classifica- 
tion scheme, the detected HI lines are type I profiles, i.e. they 
are nearly symmetric showing no evidence for absorption fea- 
tures, but they have large FWHM and extended wings. All the 
Paschen and Brackett lines detected, in fact, are broad (AV~200 
kms except for the Br 13 and the Br 15 lines for which we 
get only a lower limit to the FWHM) and slightly blueshifted 
(-5 < Vpeak < -30). In particular, the intensity profiles normal- 
ized to the continuum emission of the Pay and Pa/3 lines are 
shown in Figure |2] Both the Pay and Pa/? lines have a broad 
nearly symmetric profile peaking at 3-3.5 times the intensity of 
the continuum. The profile is centered at -5 and -11 kms re- 
spectively, and extends at least from -300 km s"' to +200 km s"' . 

The ratios between the lines of the Paschen and Brackett se- 
ries can be used to retrieve further informations about the origin 
of the HI lines. The NACO spectra have too low spectral reso- 
lution to derive the properties of the line profiles (Vpeak, AV). 
They can be used, however, to measure the equivalent widths 
and the fluxes of the Paschen and Brackett lines (see Tab. |2]). 
With respect to the ISAAC data, these spectra have the advan- 
tage of covering all the Brackett lines from the Br 10 up to the 
Br 19 in the same spectral segment, allowing us to measure and 
compare the fluxes of all the lines without introducing calibra- 
tion problems. First, we corrected the observed fluxes for the 
reddening. We used the average value of the visual extinction 



estimated by iGiovannelh et al.l ( Il995h in the circumstellar re- 
gi on of RU Lupi (Av= = 0.7 + 0.3) and the extinction law derived 
bv lRieke & Lebofskvl(ll985h . The dereddened line fluxes of the 
Paschen and the Brackett lines are summarized in Tab.|2]for the 
NACO data and in Tab.[3]for the ISAAC ones. The error aff^ect- 
ing the line fluxes and their ratios is computed considering the 
noise over the continuum and the uncertainty on the Ay value. 
Then we computed the Brackett decrement, i.e. the ratios be- 
tween the intensity of the Brackett lines coming from the upper 
level n and the first line of the series. The first Brackett line in the 
spectral range covered by our spectra is the Br 10, but this line is 
not well corrected for the atmospheric transmission (see Sect.|2]i. 
Thus we used the ratios with the Brackett 1 1 line. The Brackett 
ratios, Brn/Brll, as a function of the upper level n, are plotted 
in Fig.[3](the ratios Brl5/Brl 1 and Brl3/Brl 1 are not taken into 
consideration since for the Br 13 and Br 15 lines we have only an 
upper limit of the flux, see Sect. |2]l. In the Figure we overplot- 
ted the theoretical values of the Brackett ratios calculated for 
two diff'erent cases: (i) the solid line traces the ratios in the case 
of optically thick emission at the temperature T=10'* K (black- 
body emission) ; (ii) the four dotted lines correspond to Case B 
recombination jHummer & Storevill987[) for T^IO'^ K and four 
different values of the electron density. Note that for the consid- 
ered values of the electron density, variations of the temperature 
between 10^ and 10^ K leave the curves almost identical. The 
plot of Fig.[3]shows that the ratios between the Brackett lines are 
well fitted by case B recombination, i.e. by optically thin emis- 
sion. As we will discuss in Sect. 14.11 however, case B cannot 
reproduce the observed fluxes for typical values of the emitting 
volume. 
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Fig. 3. Ratios of the Brackett lines with respect to Brll versus 
the upper quantum number n. The solid line indicates the the- 
oretical value of the ratios expected in case of optically thick 
emission at T=10000 K. The four dotted lines are the ratios ex- 
pected from case B recombination (Hummer & Storey 1987) as- 
suming T^IOOOO K and n,=10\ lO'*, 10^ 10*' cm"^ going from 
the lower to the upper line. 



level becomes easily over-populated. The absorption feature is 
generated mainly via resonant scattering since the /} 10830 transi- 
tion is the only permitted radiative transition from its upper state 
to a lower one and the electron density is unlikel y high enough 
to cause col lisional excitation or de-excitation. lEdwards et al.l 
(I2003L |2006|) demonstrate that the shape and width of the Hel 
/il0830 line can constrain the wind launch region. They show 
that a broad and deep blue absorption, as the one detected in our 
spectra, require formation in a wind emerging radially from the 
star rather than from the disk. In fact, the broad range of veloci- 
ties covered by the absorption feature can be explained by radi- 
ally emerging wind streamlines absorbing the 1 fj.m continuum 
from the stellar surface, thus tracing the full acceleration region 
of the inner wind. On the contrary, a disk- or x- wind is con- 
fined to nearly parallel streamlines emerging at an angle to the 
disk surface and thus the continuum photons from the star will 
intercept a narrower range of velocities giving r ise to narrow ab- 
sorption feature as the ones shown in .Edwards et al. 1 (12006.) . In 
conclusion, the absorption feature in our spectra appears to in- 
dicate the presence of an inner wind radially emerging from the 
star Further insight on the geometry and the presence of differ- 
ent wind contributions to the Hel /110830 line can be retrieved 
from the spectro-astrometric analysis presented in Sect. [33] and 
discussed in Sect. 14.21 

The emission part of the line, instead, is produced by the 
accreting gas columns and, partially, by the resonant scattering 
process. The contribution from the latter depends on the electron 
density and in our case should be much lower than the loss of 
flux in the absorption part of the line. In fact, the detection of 
the absorption feature indicates that the efficiency of collisional 
de-excitation is much lower with respect to the radiative decay. 



To avoid systematic errors that can affect the lines of the 
same series, and to constrain even more tightly the emission 
conditions, we considered also the ratio between the observed 
Paschen lines, Pa/6/Pay. The measured Pa/S/Pay ratio retrieved 
from the NACO spectra is ~0. 8+0.1, and the one retrieved from 
the ISAAC ones is ~ 1+0.1. Such values cannot be produced by 
the case B recombination. In fact, for any choice of the electron 
density and the temperature between 10^ and 10'" cm"^ and 5 
10^ and 3 10"* K, respectively, case B gives Pay6/Pay >1.4. On 
the other hand, in the case of optically thick emission (black- 
body emission) the Pa/J/Pay ratio is ~0.6 at T=10"' K, and ~0.7 
at T=5 10^ K. This implies that the Paschen lines are not repro- 
duced by one of these limiting cases and that their ratio strongly 
depends from their relative optical depth. 

3.2. The Hel A1 0830 line. 

The Hel line profile obtained from medium resolution ISAAC 
spectra (Fig. |2]l shows a typical P-Cygni profile with the emis- 
sion part peaking at +70 kms"', and extending up to at least 
+430 kms"', and a broad blueshifted absorption penetrating up 
to 80% of the continuum at V — 200 km s 'and ranging from - 
360 kms"' to -90 km s"' (see Tab. [3j. The absorption feature 
clearl y indicates the presence of an inner wind (Edwards et al., 
l2003h . The efficiency of the Hel ^10830 transition in tracmg 
outflowing gas is due to the "metastability" of its lower level 
(2i-'S). Even if energetically far above the singlet ground state 
(21 eV), this level is radiatively isolated from it and collision- 
ally de-excitation is very low. This means that, once it is pop- 
ulated, via photoionisation followed by recombination and cas- 
cade and/or via collisional excitation from the ground state, this 



3.3. Spectro-astrometry 

Further to the standard spectroscopic analysis, we performed a 
spectro-astrometric analysis of the data. Following this proce- 
dure the displacement of the emission centroid of the li nes with 
respect to the source continuum is determined (e.g., iBailevl 
1998). This means that one can detect extended emission in a 
line beyond the seeing limit, down to (sub-)AU scales. On the 
other hand, the technique is sensitive only to asymmetric emis- 
sion and, since it measures a barycentre, it gives only a lower 
limit to the size of the emission region. The technique con- 
sists in accurately measuring the trace of each A-B spectrum by 
Gaussian fitting the spectrum spatial profile at each wavelength. 
In this way one recovers the wavelength dependent centroid 
(i.e. the so called "position spectrum") and gaussian FWHM 
of the spectra. The residual FWHM, obtained subtracting the 
FWHM of the continuum emission, can be sensitive to differ- 
ential size variations in the emission region. Ex amples of the 
applic a tion of spectro-astrom etry can be found in iTakami et al.l 
(l200l . lWhelan et alJ(l2005l) . 

The essential requirements for this kind of analysis are 
a good spectral and spatial sampling, and a high signal-to- 
noise ratio (SNR). Thus we applied spectro-astrometry only to 
the ISAAC spectra, in which the lines are very well resolved 
(SNR~100) and wefl sampled spectrafly (AV~200kms"', spec- 
tral sampling -14 kms"'), while the spatial sampling is mod- 
erately good (FWHM~0'.'6-a.'8, spatial sampHng: a.' 147). Since 
the trace of the star is tilted in the ISAAC spectrograph, the cen- 
troid is corrected for this distortion with a polynomial fit. The 
typical spectro-astrometric error is given by the accuracy in mea- 
suring the centroid of the spatial profile through the gaussian fit 
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Fig. 4. ISAAC spectroscopy of RU Lupi Pay (left) and Hel (right), with the sHt parallel to the jet and perpendicular. The line profiles 
normalized to the continuum emission are presented in the upper panels - A and A' . The spectro-astrometric signal is presented in 
the second panels - B and B'. The overplotted solid line is the spectro-astrometric signal due to a distorte d PSF caused by tracking 
error or unstable active optics. The signal is simulated following the procedure of 'B rannigan et al.] (1200611) (see text). The residual 
signal obtained subtracting the simulated signal to the detected one is shown in the third panels - C and C. Finally, the residual 
FWHM of the gaussian fit is presented in the lower panels - D and D'. The vertical dotted lines show the zero velocity of the line. 
The horizontal dotted lines indicate the continuum emission level and the +3cr threshold due to the signal-to-noise. 



and is equal to cr~0.5 FWHM/ ^/A^ dTakami et all |2003|) . For 
the considered lines {Npi,>10^), this accuracy turns out to be ~1 
mas. In the case of the He line, however, the signal is fainter in 
the absorption part {Nph~lQ*-lQ^)- Therefore, we adopted cr~3 
mas, which is the average value of the error across the line pro- 
file. Moreover, since the spatial sampling for the ISAAC data is 
not exceptional, we used a conservative value of the threshold of 
3cr (see Fig.HJi. We checked the positio n spectra for the presenc e 
of artifacts following the procedure of iBrannigan et al] (l2006h . 
These authors showed that the distortion of the PSF caused by 
tracking errors of the telescope or unstable active optics dur- 
ing an exposure, can induce artificial signals. These instrumen- 



tal biases can be simulated in shape and magnitude calculating 
the spectro-astrometric signal obtained f rom a modeled distorted 
PSF. We used the same procedure than in lBrannigan et al.l (l2006h 
to simulate the biases that can be produced with our instrumen- 
tal setting (cf. Tab.[T]). The shape of the artifacts produced by a 
distorted PSF depends on the slope of the line profile. In partic- 
ular, the peaks of the simulated biases appear at velocities where 
the gradient of the line profile is the largest. This explain the 
double peak of the simulated spectro-astrometric signal in the 
blueshifted part of the Hel line (see panel B and B' of Fig.|4]l. 
The magnitude of such biases, instead, depends on the amplitude 
of the distortion of the PSF with respect to the slit axis, i.e. on 
the parameters Axi, Ax2, Ayi, Ay2 in Eq. (4) of IBrannigan et al.) 
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(l2006h . The maximum biases are obtained when the distortion 
of the PSF in the slit is maximum, i.e. when we put in Eq. (4) of 
[Branniganet al. (2006) A.^:i = Ax2 = O'.'IS, Ayi = Ayi = O'.'IS, 
so that the total slit illumination displacement {Ay) is equal to the 
slit width, 0'.'3. The resulting spurious signal was then subtracted 
from the position spectra. 

The position spectra obtained for the Pa/? and Brackett lines 
do not show any shift of the centroid with respect to the contin- 
uum emission. On the other hand, in the Pay and the Hel line we 
found a spectro-astrometric signal (see Fig. |4]i. In panels A and 
A' the intensity spectrum of the two lines in the two slit positions 
is shown. The position spectrum, in panels B and B', indicate a 
shift of the centroid of the emission with respect to the contin- 
uum emission both in the parallel and in the perpendicular slit. 
The superimposed solid line is the spurious spectro-astrometric 
signal that can be induced by an elongated PSF in the slit in 
the worst possible case. In panel C and C the residual signal 
obtained subtracting the determined biases is shown. Finally in 
panels D and D' we plot the residual FWHM of the gaussian fit. 

Fig. |4] shows that there is no real signal in the Pay line. In 
fact, there is no residual in the position spectra of this line after 
the subtraction of the biases as well as in the spatial FWHM 
across the profile. In the Hel line, on the contrary, the posi- 
tion spectra cleaned by artifacts (panels C and C) show a shift 
over the signal-to-noise limit (9 mas, corresponding to ~1.3 
AU) in the blueshifted part of the line. The residual FWHM 
across the profile presents a significant signal in the same veloc- 
ity range where the spectro-astrometric signal is detected (-50 
kms~'<V<-200 kms"'), reaching differential sizes of 0'.'025. A 
broader FWHM is expected in the region where the absorption 
feature is detected, because the continuum emission is much 
fainter. It is, however, interesting to notice that: (i) the spatial 
profile can be identified very well over the background noise 
across the entire Hel spectral profile (the signal-to-noise ratio is 
S/N~100 in the spectral region covered by the absorption fea- 
ture and ~10 times higher in the emission part of the line); (ii) 
the broadening of the line is not symmetric with respect to the 
peak of the absorption feature {'Vpeak~-20Q kms"'). Instead, a 
constant or increasing FWHM in the parallel and perpendicular 
slit respectively, is measured in the region where the spectro- 
astrometric signal is detected (-200 kms"'<V<-100 kms"'). 
This suggests that multiple spatial components are contributing 
to the line profile in such region of the "P Cygni" absorption fea- 
ture. The interpretation of the detected spectro-astrometric signal 
is discussed in Sect. 14.21 

4. Discussion 

As illustrated in the Introduction the excitation mechanism of 
the HI lines is an actively debated topic. Up to a few years ago, 
both the observations and the models focused on the Balmer 
lines, and in particular on the Ha line. The large amount of op- 
tical observations evidenced the limits of both the accretion and 
the wind models in reproducing the Balmer lines and motivated 
the develo pment of hybrid mode ls with more complex geome- 
tries (e.g., iKurosawa et al.Ll2006l) . On the other hand, the capac- 
ity of these models to reproduce the profiles and fluxes of the 
HI infrared lines was poorly tested. Interestingly, the observa- 
tions of Paschen and Bracket t lines, such as the o nes presented 
injFolha & Emerson (2001) jNisini et al.1 (|2004|) . IWhelan et af] 
( 12004)) . lEdwards et al.l HlOOd ), have shown that the magneto- 
spheric accretion models cannot fully reproduce the character- 
istics of the observed profiles, in favour again of an hybrid 
wind/accretion model. 



As we have shown in the previous section, the Paschen and 
the Brackett lines in our spectra present blue-shifted broad sym- 
metric profiles (y peak ~ -17 kms"' and AV ~ 200 kms"'). These 
are difficult to reproduce if one adopts separately a wind model 
or, alternatively, a magnetospheric accretion model. Wind mod- 
els generally predict red-shifted peaks and strongly asymmetric 
profiles with P-Cygni absorption features, as the one detected in 
the Hel line (Calyet et al., 1992). Hartmann et al. (1990) were 
able to reproduce more symmetric Bry and Bro- profiles, with- 
out the P-Cygni absorption feature with a model of magnetically 
driven wind, but their peaks are still red-shifted. lEdwards et al.l 
(Il987h showed that an optically thick circumstellar disk may 
cause the occultation of the red-shifted part of the profile thus 
explaining blue-shifted peaks. In any case, asymmetric profiles 
are expected. 

On the other hand , magnetospheric accretion models 
(iMuzerolle et al.L 11998a ') predict blue-shifted peaks, but asym- 
metric and narrow profile s (AV ~ 100 kms"') oft en showing in- 
verse P-Cygni absorption. lMuzerolle et al.l(l200ll) retrieved from 
their accretion model more symmetric and broad line profiles, 
taking into account line damping due to different broadening 
mechanisms, but this effect is really important only for Ha and 
Balmer lines, while it is negligible for Paschen and Brackett 
lines. 

We will now examine how the above scenarios fit in the case 
of our RU Lupi spectra. 

4. 1 . Accretion vs wind 

As illustrated in Sect. 13. li the case B recombination well repro- 
duces the Brackett decrement, but not the Pa^S/Pay ratio, whose 
value implies that these two lines are optically thick. Moreover, 
case B recombination assumes (i) a constant density, which is 
not applicable to the circumstellar region of T Tauri stars, and 
(ii) optically thin lines. To reproduce the observed fluxes, how- 
ever, the latter hypothesis would imply an emitting volume cor- 
responding to a spherical envelope extending up to 40 AU from 
the source, which is unrealistic according to the models of cir- 
cumstellar accretion/ejection structure (e.g., Nisini et al., 2004|). 
To investigate the origin of the HI lines one has to assume a 
variable density profile, and the optical depth and the flux of the 
Paschen and Brackett lines will be retrieved integrating along the 
line of si ght. 

Nisini et al.l (|2004|) analysed the Brackett lines in the Class I 
source HH 100 IR, finding line profiles very similar to the ones 
of our Paschen and Brackett lines. Interestingly, they found that 
the Brackett decrement in HH 100 IR is neither reproduced by 
the Case B recombination curves (i.e. optically thin emission), 
nor by the blackbody emission (i.e. completely thick emission). 
The authors suggested that these lines are emitted in an expand- 
ing envelope of ionised Hydrogen, and thus different lines are 
tracing different layers of the envelope according to their optical 
depth. Following this idea they were able to reproduce the ra- 
tios between the optically thick Brackett lines through a simple 
model of spherical wind. 

We used the same approach to reproduce the Brn/Brl 1 ratios 
(Fig. O, the Pay6/Pay ratio and the observed line fluxes both in 
the accretion and in the wind hypothesis adopting a toy model 
consisting in a spherical envelope of ionised Hydrogen, extend- 
ing from an inner radius r, to an outer radius rout, where LTE con- 
ditions are assumed. If the envelope is expanding (wind hypoth- 
esis) or infalling (accretion hypothesis) the radiation transfer can 
be treated in the Sob olev approximati on followin g the formalism 
described in lCastod ([T97Q) and iNisini etal.1 ([T995ft . The gas can 
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be only partially ionised, and, therefore, the envelope total mass 
loss or accretion rate is Mwmdjaccr - MioiJx,,, where is the 
ionisation fraction. The electron density in the envelope is given 
by the continuity equation, - MKi„diaccrXe I ^mnHr^Vir). 

Using this model both the wind and the accretion hypothesis 
can be tested just changing the velocity law V{r). In fact, this 
will imply a different density profile and thus a difference in the 
relative optical depths of the lines in the two cases. 

In the wind hypothesis the HI envelope is expanding and the 
gas velocity follows a general law of the type V{r) - Vo+{V„,a.x- 
Vo)(l -('"i/'")"), i-e. the gas is accelerated from an initial velocity 
at the base of the wind, Vo, at the maximu m velocity Vmax at a 
distance that depends on the parameter a (iNisini et all l2004t) . 
Vmax is the maximum velocity of the expanding gas. The FWHM 
of the observed lines is determined by the velocity of the gas 
particles, thus we retrieve from the observed profiles V,„„v~200 
kms"'. 

The simulated emission is tuned by a proper choice of the 
input parameters. For typical values of the size of the envelope 
(r,- = 1 Rt, rout = 3 - 10 R,) and of the electron temper ature (T 
= 10^ K), and taking a-4 and V'o=30 kms 'as in Nisi ni et aTl 
we are able to reproduce the Br«/Brll ratios and the 
absolute fluxes of the Paschen and almost all the Brackett lines 
(within +4cr) for M,o„ ~ 2 10"** Mq yr"'. On the other hand, in 
this way we are not able to reproduce the observed PayS/Pay ra- 
tio (Pa/S/Pay~1.3-1.4). Considering a slower acceleration of the 
wind (a-l-3) and/or a lower temperature (T down to 5 10^ K) 
the Pay6/Pay ratio decreases to ~1-1.3 and the absolute fluxes 
are in agreement with the observed ones but the Brackett ra- 
tios are not reproduced because the lines optical depths are too 
large. Taking lower values of the mass loss rate (M,o„ < 10"** Mq 
yr"'), that is the fundamental parameter regulating the emission 
in the HI lines, we obtain a good fit of the Brackett decrement 
for rout within 3 - 10/?* and values of the exponent in the ac- 
celeration law of the wind between 1 and 4. The absolute fluxes, 
however, are too low in this case (3-4 times smaller than the 
observed ones) and the Pay6/Pay ratio is >1.4. To obtain the re- 
quired Pa/3/Pay one should assume a value of M,o„ higher than 
2 10"** M0 yr"'. For such a massive wind, however, the Brackett 
decrement and line fluxes would no longer be fitted by the model 
because their optical depth would be too large. 

Since we could not find a choice of the parameters for the 
wind model reproducing at the same time the Brackett decre- 
ment and the observed Pa/?/Pay ratio, we tentatively tried to use 
the model to simulate the emission from a spherical accreting 
envelope, by reversing the sign of the gas speed. A spherical 
model is not the best approximation of the accretion process, 
that, according to the most recent models, proceeds along the 
magnetic field fines, in a n axisymmetric non-sp herical geometry 
dHartmann et all 119941: iMuzerolle et all Il998a ). Nevertheless, 
we have adopted the spherical model as a first approximation, 
in order to check if the observed ratios are better reproduced by 
accretion rather than by an outflow. 

If we assume that the HI lines are mainly excited in the ac- 
cretion columns we can derive an estimate of the mass accretion 
rate from the luminosity of the Pa/3 line. To this aim, we use the 
empirical relation between L(Paj3) and the accretion lumi nosity 
derived bv IMuzerolle et all (Il998d) and lNatta et all (|2002|) : 

log L„cc/Le = 1.36 \ogL{Pal3)/Ls + 4. (1) 

The mass accretion rate is then computed from Lace as M„oc - 
LaecRJi GM,). Using for RU Lupi M, ^ 0.8 M p and R, = 
1.6 Rq dLamzin et all 119961: iHerczeg et all l2005b . we obtain 
Mflcc=2 10"^ M0 yr"'. This value is in agreement with previous 



estimates of iLamzin et al.l (Il996h . that tuned the value of Mace 
in order to reproduce through their model the observed contin- 
uum energy distribution (see, however, Herczeg et al. 2005, for 
an alternative derivation of Mace using the accretion luminosity 
estimated from the UV excess). 

In the case of an accreting envelope the typical velocity law 
is V{r) = Vniaxf '^'^- We fixed the size of the envelope and the 
temperature (r, = Rt, rout = 3/?», T^IO"* K), and examined the 
results obtained varying the ionisation fraction x^., and thus the 
electron density of the accreting columns for Mace-"^ 10"^ M© 
yr"'. If the gas is fully ionised (Xf=l) we obtain Paj0/Pay~O.7, 
but the Brackett decrement is not reproduced, because the opti- 
cal depth of these lines is too large, and the absolute fluxes of all 
the lines are higher than the observed one by one order of mag- 
nitude. Indeed, estimates of the ionisation fraction for different 
values of the temperature and the mass loss rate show that the 
gas is only partially ionised in the circumstellar envelope (e.g., 
iNatta et al.l 11988) . Assuming Xe-OA both the Brackett decre- 
ment and the ratio between the Paschen lines are well repro- 
duced, but the absolute fluxes are still too high by a factor 3-4. 
To obtain absolute fluxes in agreement with the observed ones 
we have to assume low ionisation fractions, down to Xe=0.15. 
In this case, however, we are not able to reproduce the Pa/3/Pay 
ratio. We note that in a non-spherical geometry the lines become 
optically thick for lower mass accretion rates, and thus both a 
"thick" Pa/S/Pay ratio and lower absolute fluxes can be obtained 
at the same time. In such a model, however, the profiles of the 
HI lines are predicted to be much narrower than the observed 
ones, i.e. FWHMs of the order of 100 km s"' would be obtained, 
in contrast with our observations. 

In conclusion, our analysis shows that neither a spherically 
symmetric wind nor spherically symmetric accretion can repro- 
duce the Brackett and Paschen line ratios and fluxes. There are a 
few recent observational results suggesting that both the emis- 
sion from the accretion colum ns and from a w ind may con- 
tribute to the line profiles. Takamfet alj(l2001h and lWhelan et al.l 
(|2004|) . for example, found that the gas from a wind contributes 
to the emission in the wings of their He and Pa/3 lines, thus 
producing line profiles broader than the ones predicted by mag- 
netospheric accretion models. In our case, however, the observed 
H lines do not show the core/wings structure which is expected 
when both the accreting material and a wind are contributing to 
the emission. 

To investigate further the possibility of a co-existence of 
emission from both inflowing and outflowing gas in our H 
lines we analysed critically the spectro-astrometric results as de- 
scribed below. 

4.2. Insights from spectro-astrometry 

In the last two Sections we showed that both the accretion model 
and the wind model fail in reproducing all the characteristic of 
the observed NIR HI lines. A powerful technique to investigate 
the origin of the HI emission and separate the accretion and the 
win d contributions to the line profiles is spectro-ast rometry (see, 
e.g.. lTakami et al. L l2001l:IWhelan et alil2004l2005h . In fact, out- 
flowing gas presents an extended emission testified by a de- 
tectable positional shift, while accreting gas remain concentrated 

on the source location. 

Although lTakami et al.l (I2OOII) detected extended emission in 
the wings of the Ha line up to ~3 AU from RU Lupi, we found 
that there is no spectro-astrometric signal in the HI NIR lines 
(see Fig.Hl). The lack of signal indicates that: (i) either the emis- 
sion comes from a compact region (smaller than ~1.3 AU); (ii) 
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or, if there is a further emission from an extended region, this 
is symmetric, thus not detectable through spectro-astrometry, or 
too faint with respect to the emission from the star and the ac- 
cretion columns (see below). 

On the contrary, we do detect after subtraction of possible 
biases a substantial spectro-astrometric signal in the Hel /110830 
line (Fig.m panels C, C). The shift found in the Hel line is inter- 
preted as an emission originating in the base of the jet, that par- 
tially fills the absorption feature of the "P Cygni" profile due to 
the inner wind. The inner stellar wind acts like a natural corono- 
graph hiding the star along the blueshifted absorption, increas- 
ing the contrast of emission from outer regions. This could ex- 
plain why we detect a positional shift in the Hel line and not in 
the NIR HI lines. A similar extended emission superimposed to 
the blueshifted absorption feature of the Hel line was found by 
[Takami et al. (2002) in DG Tauri. Following the same reason- 
ing of iTakami et al.l (|2002|) we suggest that there are three con- 
tributes to the profile of the Hel line: (i) red-shifted emission pre- 
sumably due to accreting gas; (ii) blue-shifted absorption due to 
the presence of an inner spherical wind emerging from the stel- 
lar surface; (iii) blue-shifted extended emission between V — 50 
kms"' and V~-200kms ' superimposed to the absorption fea- 
ture and coming from a jet at a larger distance from the source. 
This emission is supposed t o come from t he ba se of the colli- 
mated micro-jet detected by iTakami et al.l (120011) in the forbid- 
den [S ii] and [Oi] lines and in the Ha line. The velocity range 
covered by the Hel extended component is the same as in the Ha 
line. The Hel spectro-astrometric signal extends up to ~20 mas 
from the source (~3 AU). However, as explained in Sect. 13.31 
this is only a lower limit to the true angular scale of the emission 
region which could be more extended. 

Combining in a bi-dimensional plane the shifts detected 
along the parallel and the perpendicular slits we find that the 
extended em ission is not completely aligned with the jet of 
ITakami et al I (l200lb . the misalignment being not greater than 
30°. This may be justified by the fact that: (i) the jet may be not 
well collimated in the He line; (ii) our observations were made 
4 years later and a 200 kms~' s tream moves 1 AU in a week. 
Furthermore, Taka mi et al.l (12001 ) do observe variability in their 
spectro-astrometric data taken several years apart. Another dis- 
crepan cy with respect to the result obtained by I Takami et al.l 
( 1200 lb in the Ha line is that we did not detect extended emission 
in the red-shifted part of the Hel line. This can be due, again, to 
the fact that the jet emission in the HI and Hel lines is fainter 
than in the optical Ha line and it can be evidenced by spectro- 
astrometry only over the blueshifted absorption, which hides the 
emission from the star and increases the contrast of outer emis- 
sion. 



5. Conclusions 

In this paper we analysed infrared spectra and images of the T 
Tauri star RU Lupi, in order to find observational constraints for 
the physics of the circumstellar region. This source was previ- 
ously observed in the optical wavelength range by other authors, 
that reported about a strong accretion activity and the presence 
of a mic ro-jet detected in th e optical forbidden lines and in the 
Ha line (iTakami et al.Ll200ll) . We have observed for the first time 
RU Lupi in the NIR with high angular and spectral resolution. In 
this paper we report about these observations and, in particular, 
we investigate the origin of the permitted H and He NIR lines 
through the analysis of the line profiles and fluxes, and using 
the spectro-astrometry technique which allows us to constrain 



the emission region of each velocity component of the detected 
lines. 

The problem of the origin of HI lines was widely de- 
bated in the last decades. Many different models were pro- 
posed to reproduce their profiles and fluxes: w ind models (e.g., 
Hartmann et al., 1990), accretion m odels (e.g., iHartmann et al.L 
1994; Muzero Ue et al.(.il998all200 1) and also hybrid models ac- 
counting for the emission from both the accretion columns and 
the wind (e.g., iKurosawa et al.L l2006l) . These analyses, how- 
ever, focused mainly on Balmer lines observed in the optical 
range. The Paschen and Brackett lines detected in our NIR spec- 
tra of RU Lupi present broad, slightly blueshifted and nearly 
symmetric profiles, which are difficult to reproduce with either a 
win d or an accretion model. This was already note d , in general, 
bviFolha&EmersonI (1200 lb . IWhelan et al] (|2004|) . iNisini et al.l 
(12001 7 We tentatively used a toy model of a spherical envelope 
of partially io nised Hydrogen w ith a wind or an accretion ve- 
locity profile dNisini et al.L 120041) to constrain the HI line ratios 
and fluxes. This analysis showed that neither a spherical wind 
nor spherical accretion can reproduce the Brackett and Paschen 
decrement and fluxes. 

Our spectro-astrometric analysis did not highlight any ex- 
tended emission in the Paschen and Brackett lines, suggesting 
that the region emitting the HI NIR lines is very compact (<1.3 
AU) and/or symmetric. On the other hand, the Hel /110830 line 
blueshifted absorption feature clearly indicates the presence of 
an inner wind. The wide range of velocities covered by the ab- 
sorption feature favours the geometry of a spherical wind emerg- 
ing from the st ellar surface ra t her th an a disk wind as found in 
many TTS by lEdwards et al.l (12006*). In addition, the spectro- 
astrometric analysis highlighted the presence of an emission su- 
perimposed to the absorption feature that extends up to at least 
~3 AU from the source and covers the sam e velocity rang e 
of the Ha extended emission found by Tak ami et al.l ( 1200 lb . 
We suggest that this emission comes fro m the blue lob e of th e 
micro-jet detected in the optical lines by iTakami et al.l (1200 lb . 
Interestingly, the Hel showed to be more sensitive than the HI 
NIR lines to faint extended emission because the absorption fea- 
ture increases the contrast between the emission from the ex- 
tended region and the continuum from the source. This confirms 
the potential of the Hel /il0830 line to investigate the complex 
geometry of the inner part of the circumstellar region of CTTS, 
where both the accretion and the ejection processes take place 
(see also Edwards et al., 2006). On the base of these results, 
future high angular resolution observations of the permitted H 
and He lines in CTTS can give useful hints to the understand- 
ing of the accretion/ejection mechanisms and the development 
of theoretical models accounting for the emission from both the 
outflowing and the accreting gas. 
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